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Abstract 

o ; 

^ Models of PT symmetric quantum mechanics provide examples of biorthogonal quantum systems. 

The latter incorporporate all the structure of PT symmetric models, and allow for generalizations, espe- 
fT^ ' cially in situations where the PT construction of the dual space fails. The formalism is illustrated by a 

^ , few exact results for models of the form H = {p + v)'^ + X]fc>o ^-'^P {i^^)- In some non-trivial cases, 

■ equivalent hermitian theories are obtained and shown to be very simple: They are just free (chiral) 

^ _ ^ particles. Field theory extensions are briefly considered. 

I> 

o 
o 



^ • *curtright@physics. miami.edu & mezincescu@physics.miami.edu 

+-> 

c3 ■ 1 Introduction 

O^' There has been some recent theoretical interest in non- hermitian Schrodinger equations, especially in the 

guise of "PT symmetric theories" Therefore, we have considered a few simple soluble examples and 

• ■ explored them in some detail. We believe exact solvabihty permits the underlying structure to be appreciated 

r> I more completely. To this end we have found Neumann's work |70| to be remarkably prescient. Although 

■ in physics there has long been an interest in complex Hamiltonians for phenomenological purposes 41,, 

I4()| . quite generally mathematical interest in non-self-adjoint spectral problems pre-dates the invention of 
quantum mechanics jl2| . See j37 l for a brief survey of the early mathematical history. 

Some previous papers have touched on the relevance of biorthogonal systems for PT symmetric models 
|59[ I60L IHTl 1671 188| . Here we wish to stress the importance of such systems, and their generality, in the 
context of elementary soluble cases. We consider models whose Hamiltonians are of the form |53[ 1731 1421 

EnmsiiniEiiiizi 

H — {p + + ^ exp {ikx) (1) 

where the imaginary exponents are of the same sign, and where v and /ifc have arbitrary values. Although 
much is known about the mathematical structure of these models, they tend to be overlooked in the physics 
literature on PT symmetric theories. We seek to remedy this by offering a thorough discussion of their 
properties, especially in the case of periodic solutions, for which we believe we make several new observations. 
We analyze at considerable length the cases of one and two exponentials, and then we indicate how the results 
extend to the general case. We believe it is almost impossible to over-emphasize the merits of going through 
such exactly solvable models in complete detail - "the simpler the better" if the important structural features 
are retained. We end with a brief look at field theoretic extensions. 
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2 Biorthogonal systems 



A sequence of elements {tpj} and linear functionals {A^^} is said to be biorthogonal (or more precisely, 
hiorthonormal) when 

Afe ii^j) = (2) 

See the classic texts by Goursat (021 Vol. Ill §573 and §594), Banach (jH] Chapter 7), Morse and Feshbach 
(inni §7.5 and §8.3), and Gohberg and Krein (gj Chapter 6), as well as [SI H EHl HI EDI ■ The notion 
of biorthogonality is especially useful when dealing with elements that are analytic functions defined on a 
contour in the complex plane. In this case, the linear functionals are usually realized as integrals along the 
contour involving a set of "dual functions" Xk = i''Pk)^^^\, so 



^k (ipj) ^ J Xk (z) Tpj (z) dz = Sjk (3) 

However, it is important to stress that a priori there need not be a simple point-by-point relation between 
the functions and their duals (unlike the standard situation in quantum mechanics). In general, Xfe {z) is 
not obtained just by conjugation, and differs substantially from ipj (z), even when the contour of integration 
is a real line segment. Moreover, the normalization and phase of tjjj can be changed, arbitrarily, if a 
compensating change is made in the normalization and phase of Xj- That is to say, if the system {ipjjXk} 
satisfies ©, then so does the system {Zjipj,Xk/Zk} (no sum on j, k) for all non- vanishing, finite choices of 
the Zs. 

For any function in the span of {V'j}: we have 

tp{z)^y^ Cnlpn (z) , C„ = / Xn (z) Ip (z) dz (4) 



assuming the requisite convergence of the sum. For such a -0 [z) we may define, at least formally, an 
associated dual function by 

V'dual (z) = ^ clXk (2) (5) 
k 

While the convergence of the sum for ^/^duai {z) is not necessarily guaranteed by the convergence of that for 
ip (z), nevertheless, if it is permissible to interchange and J dz, then we can make use of (gj and define 
a norm as 

II VII' = f ipdu.iiz)iPiz) dz = ^|c„|'>0 (6) 

•' n 

If is finite, this allows us to define a bona fide probability distribution on n corresponding to the 

fimction tp. The probability of n is just |c„| / \\Tp\\ . However, it is important to stress that the bilinear 
ipdimi (z) ip (z) might not provide a bona fide probability distribution in the variable z, since i/'duai (z) ip (z) 
might be complex. Moreover, J V'duai (z) z %p (z) dz might be complex for one or both of two reasons: z 
might be on a complex contour, and ■0duai (z) is not necessarily ip (z), in general. (Perhaps it is possible to 
think of V'duai (z) ip (z) as a "complex probability" density but we will not pursue that idea here.) 

When n labels the eigenvalue of an operator L, such that Lipn (z) = Xnipn (z), we may also compute 
averages of functions / (L) using a procedure analogous to that of conventional quantum mechanics, even if 
"Vduai (z) ^ ijj (z). These statistical averages are defined to be 

if)i, [ ^dual (Z) / (L) 4, (Z) dz^Y. / (^») / ll^ll' (7) 

IIV'II J 

again assuming the requisite convergences. Whether these averages are real depends on both / and {A„}. 
Also, for any given L, the action of the operator on V'duai (z) is not necessarily obvious from its action on ip. 
Similar statements apply if we interchange the roles of ip and f/'duai- 

For any biorthogonal system, it is always possible, at least formally, to express the norms © for functions 
in the span of {ipj} as integrals of bilinears in ip and ip* , where the construction is bilocal, in the general 
situation. If we build a kernel from the dual basis functions as follows 

K {w,z) = YXn{w) Xn{z) (8) 
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once again assuming the requisite convergences, we then have 



V'duai (z) = / V' (w) K {w, z) dw (9) 



IIVjII = J J ^P{w) K {w,z) V(^) dwdz (10) 

where the contour for W is the complex conjugate of that for z. Similar results follow, at least formally, 
from constructing a dual bilocal kernel and integrating this with the associated dual functions. 

J {W,z) = ^lljn{w) i>niz) (11) 



■0 (w) = J (W, z) V'duai (z) dz (12) 



IIV'II^ = yy V'duai (w) J {W,z) V^dual (z) dwdz (13) 

Also, note that the convergences of the sums for J and K can be affected by making compensating changes 
in the individual normalizations of the elements in the biorthogonal system: {tpj, Xk} ~^ {^j''Pjj Xk/Zk}- By 
suitable choice of the Zs, the sum for one of the kernels can be made to converge uniformly, and rapidly. In 
the process, the sum for the other kernel may well become a divergent (perhaps asymptotic) series. 

Formal expressions for other kernels follow immediately from acting with various operators on the {ipj} 
or {xfc} functions, or their conjugates, under the sums. For instance, returning to the previous eigenvalue 
example, 

Jf{L) (w, z) 

EE ^ Vn (W) f (L) V„ (Z) = 5] V-n (w) f (A„) V« (z) (14) 
n n 

which allows writing the average |7J) as 

1 



(/)^ = / / V'duai (w) Jfi^L) {w, Z) V'duai {z) dwdz (15) 



In the Jf[L) {w, z) form, hermitian kernels leading to real averages are easily recognized when / (A„) — / (A„) 
for all n. Similar statements hold for kernels -f^£((Lj„_,i) built of bilincars in the dual functions, when the 
Xn {z) are eigenfunctions of operators Lduai- 

Examples For an elementary example, let {V'fc} = {^z^ \ fc = 0, 1, 2, • • • }, i.e. the usual basis for analytic 
functions near the origin, and take the contour to be any simple closed curve surrounding z = 0. Then 
Xk {z) = z~^~^ I^Txi, for fc = 0, 1, 2, • • • . If the contour is a circle of fixed radius i?, with z = i?exp {—iO), 
and if for convenience we incorporate various factors into a measure, dfi — —dz/2Triz — d6 /2tt, then the 
previous sequence {^'^j becomes the familiar set of left-moving eigenfunctions of pg — —id/dO, {V'fc [9)} = 
{R'' exp {—ikO) I A: = 0, 1,---}, and the dual functions become (upon removal of 1/z) the right-movers, 

{Xk{0)}^{'-^^ I A: = 0,l,---},so that 



Xk{9)i^j{0)d0 = Skj (16) 



But even in this familiar example, as defined, Xfc (0) ^ V'fc {0) unless i? = 1. The bilocal norm kernels are 
easily constructed for this elementary case: J (0, 0) = SS^Lo ^xp in{Q — Q) = 1/ (l — i?,^ exp i (^6 — 6)); 
K(9,0) = ^'^^Q R~^" exp in {^9 — 6) = i?^/ (i?^ — exp i (6* — ^) . As functions of R, the complementary 
convergence properties of J and K are manifest. 

A more interesting example was discovered by Neumann |7U) in his work on Bessel functions (see §9.1-17 
in |82|). This example is relevant to our discussion of FT invariant Hamiltonians to follow. For all analytic 
Bessel functions of non-negative integer index 



fc=0 ^ ^ 
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there are corresponding associated Neumann polynomials {An} in powers of 1/z that are dual to {Jn} on 
any contour enclosing the origin. These are given by 

Ao{z) = l, A,{z)^- An>2{z)^ni-\ ^ (-) (18) 

For convenience we have modified the usual form of the associated polynomials, 0„ (as given in |70[ I82L IT] 1 . 
and defined An (z) — Sn z On (z) where eq = 1 and Sn — 2 for n ^ 0. Upon integrating counterclockwise 
once around z = 0, 

_L A, (z) Jk (z) = (19) 

so that {Jn {z)} and {An (z)} constitute a biorthogonal system on the complex plane in the same sense as 
{z"} and {z~"}, but with the very significant difference that An (z) ^ J„ (z) even when z is on the unit 
circle. 

There are invertible maps between these two biorthogonal systems, {z^ , z"*^} < — > {Jj (z) , Ak (z)}, about 
which we will say more later. The existence of such maps can be inferred from l)17|l and H18|l as well as the 
reversions of those series to obtain all powers of z. 

l^Joiz) + 2f2j2niz) , ^/^^ = 2^- ^^±Mi|±^^— ^ Jfe+2„ (^) where fc = 1,2,3,- •• (20) 

n— 1 n—0 

1 _ 1 V {-y-'A^^jz) I _ I ^ {-Y-'A,,.^iz) 

^ (j-fc)!(j + fc)! ' z2j-i 22i-i (j - fc)! (j - 1 + fc)! ^ 

Some bilocal kernels in closed form Neumann found that the Cauchy kernel is given by (|H2Ij §9-1) 

^ - £ A„ (w) Jn (z) (22) 

When rewritten as 



w — z w 

n=0 



ld(w;,z) = ^^ = ^ J„(«;)A„(z) (23) 

n=0 

this kernel acts as the identity when integrated on the space of Bessel functions (or alternatively on Neumann 
polynomials) provided the contours of integration are appropriately chosen to handle the singularity at w = z. 

Jk iw)^^<f — Id (w, z) Jk (z) , Ak iz) = ^<f — Ak (w) Id (w, z) (24) 

^T^T'J\z\>\w\ 2 27riJ|„|<|3| w 

Acting with Bessel's operator on this kernel leads in turn to an example of a bilocal eigenoperator kernel. 

^'uP'z wz nP'z 



H (w, z) = w^-— + w— + Id (w, z) 



dw^ dw J (z-w) (z-w) z-w 

OO 

^ nV„ (u.) A„ (z) (25) 



1/" dz , , , ,,,, . o . , ^ ^ I 



k'Jk{w) = —i> —H{w,z)Jk{z), k'Ak{z) = —i — Ak{w)H{w,z) (26) 
27ri7|2|>|u,| z 27riJ|^l<l^l w 

Neumann also showed that ( 82 , §11.1-2) 

OO 

Jo (^Vw^ + Z^ — 2WZ cos 0^ — ^ Jn {w) Jn (z) COS {n(j)) (27) 
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where J_„ (z) = (—1)" Jn {z) and Jq (s) — YlV=o \k^y (^) ' ^''-'^^ (|27|l we obtain an explicit and closed 
form for the dual kernel used in the construction of the bilocal norm for this biorthogonal system, as in . 

J {w, z) = ^ J„ (w) Jn (z) = - (Jo {w-z) + Jo (w) Jo {z)) (28) 

n=0 

With the normalizations chosen, the sum for the other kernel K (w, z) — X]^o (^) (-^) asymptotic 
series. However, as discussed earlier in a general context, the J [w, z) and A" (w, z) series convergence prop- 
erties can be reversed by making compensating changes in the normalizations: {Jj,v4fc} — > {ZjJj^Ak/Zk\. 
We will give other examples of kernels, in closed form, when we discuss the corresponding quantum systems 
below. 



Generalizations The {J„ (z) ,An (z)} biorthogonal system has been generalized by Gegenbauer to other 
situations involving Bessel functions with non-integral indices and their associated polynomials (see §9.2 in 
|82j\ This is relevant to our discussion, given below, of some simple magnetic field effects. Extensions 
to confluent hypergeometric functions iFi have been proposed by Erdelyi |88j . These are relevant to 
H = {p + 1/) + fii exp (ix) + ^2 exp {2ix) for arbitrary fii and /i2, as also discussed below. We will further 
generalize the polynomials to all systems governed by Hamiltonians of the form (Q. But, before fleshing 
out the biorthogonal structure in the various soluble quantum mechanical models of interest to us, we give 
a survey of the classical dynamics for the first such model. 



3 Classical dynamics for V [x) = exp {2ix) 

Consider the Hamiltonian 

H = p^ + m^e^" (29) 

where m is real. This is a complex form of the well-studied Liouville dynamics EHI ■ The coefficient rn? 
of the exponential has to be real for the Hamiltonian to be PT symmetric. That is to say, H is invariant 
under p ^ p, x ~x, and complex conjugation. Or at least, that is the case when we restrict to real x and 
real p. These transformations are inspired by the conventional PT transformations of quantum mechanics, 
for wave functions ip (x, t) that depend on real x and t. The phase of can be changed to any arbitrary 
value just by making a real shift in x, but only shifts Ax — mT/2 for integer n will preserve PT invariance. 
Still, we may say all values of m, even if complex, give theories translationally equivalent to PT symmetric 
models. 

While this has to be the simplest periodic complex potential example one can imagine, nevertheless, it is 
not obvious a priori how or if the dynamics is consistent, nor is it obvious what connections exist between 
a classical system described by and any quantum deformation of that dynamics. A goal of this paper 
is to examine this consistency in some detail, and to discuss these connections. 

While simple, this is still an interesting dynamical system. However, there is some ambiguity about what 
H represents when p and x are complex. Following Bender et al. j6j , we choose to continue the model into 
the complex plane simply by complexification of the usual form of Hamilton's equations of motion, without 
complex conjugation of any of the variables appearing in those equations. Thus we take (|29|l and 

^^^^2p(t), 'm^JH^ -2^rn^e-W (30) 
dt op dt ax 

to hold for all complex x and p. With this definition of the continuation, the complex energy is conserved 
dH/dt = 0. So under this complexification procedure there are in fact two conserved real quantities, HeH 
and Im H. 

Complex energy conservation allows reduction to a single first order equation 

= 2p (t) = ±2\/i;-m2e2«(*) (31) 
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with integration in terms of elementary functions. This yields the classical motion 



.(t) = lln 



^cosh^ f ^2i\fEt + arctanh — ^ 



(32) 



Sign flips on the RHSs of H31|l and (|32|l occur at complex turning points, as given by Re V E — iji'2Q2ix{t) _ 
= Im \J E — m^e^"(*'. If no Re \J E — m^e^*^*^*) = points are encountered, upon taking the upper signs 
in H31|l and H32|l , the motion is a continuous but varying progression towards positive Re x with periodic 
oscillations in Imx. We plot a representative > example {E = 1/4, a; (0) = 1, Rep(O) > 0) showing 
real and imaginary parts of both x and p as parametric functions of t. 




(Rea;,Ima;) in blue and (Rep, Imp) in red, plotted parametrically for E = 1/4 



In this and all other plots to follow, we have set m — 1. Heavier/lighter curves indicate forward/backward 
evolution in t. Initial points are indicated by small boxes in the Figure. 

The momentum trajectory in the Figure is as simple as it looks: A circle. For all real E > and any 
initial x (0) the momentum p (t) evolves around a complex circle centered on the real p axis. This can be 
shown by inserting H32() into H31() and manipulating the result to obtain 

p{t)^±VEi^ ^ (33) 

where we have defined the complex amplitude 



a^-^(^l + ^l-^e-(0)j e--(") (34) 

Hence a depends on the initial data but not on t. In this form, a and p have transparent high energy limits. 
From (|33|l it is straightforward to establish by direct calculation for real E > that p (t) evolves around 

a circle in the complex plane, of radius j^^^\/E, whose center is on the real axis at zfc VE- It also 
follows from H33|) that the momentum's angular velocity around the circle is not constant, but varies. In 
addition, for real E' > a turning point is encountered if and only if \a\ = 1, implying the p circle has infinite 
radius and degenerates into a vertical straight line in the complex plane, as described by purely imaginary 
momentum trajectories p {t) = —i^/E tan (^2^/Et =F ^ arga^ . Correspondingly, the x trajectory for \a\ = 1 

is along a vertical line in the complex plane, with fixed Hex. For example, this can occur when x (0) = Ntt 
and E = m^, in which cases the x trajectory "bounces off' the real axis at t = 0. 

The E — motion is somewhat simpler, but still has interesting structure, and is given as follows. 
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Again, the momentum evolves along a circle of finite radius in the complex plane, approaching the origin as 
a limit, except in the special situations x (0) = TVtt, which correspond to circles of infinite radius. In these 
special cases, the particle stays at T{,exE=o = Nit for all but moves in the imaginary direction, with 

ImxB^o (t)lRc.,.o=w. - \ In (l + ^mh^ ± 4mi {-f) 

= ln(l±2mt) (36) 

and with ± signs depending on N as well as the choice of root for pe=q- So with the ± as exhibited in 1)36(1 . 
xe=o goes to a; = +ioo as t — s- ±cx) and goes to a; = —ioo as t — s- 1^-^- In fact, no matter what x (0) is, 
the E = solution always has an imaginary part that behaves as Imx^j^o {t) t^^^^±oo (^^) ^'^'^ least one 
choice of sign. 

We plot parametrically some E — examples, for various initial x (0), showing real and imaginary parts 
of X (t) and p{t) for forward/backward evolution. The initial positions are x (0) = 7r/2 (blue), 2.5 (green), 
2.8 (black), and tt (red), as indicated by small boxes in the Figure on the left. Corresponding initial values 
for p are also indicated by small boxes in the Figure on the right. 
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(Rex, Ima;) for E = and various initial positions (Rep, Imp) for E = and various initial momenta 



Canonical transformations Many of the classical properties are more easily understood upon taking 
note of the following fact. There exist generating functions for canonical transformations from Liouville to 
free particle dynamics |17j . with free variables 9 and pe. However, in the present context, both free particle 
and Liouville dynamics must be complexified, in general. One such transformation is given by 

d d 
F ~ me" sinff , p = —F — ime^^ sin 9 , pg = —-^F = -me" cos 9 (37) 
ox o9 

The free particle's momentum is conserved, pg = ±\^, since under the transformation 

H^p^ + m^e^" = p2 = Hi,,, (38) 

For real E > 0, pg is also real. Although this would in general require complex x, p, and 9, as follows from 
the transformation. Separating out the {x,p) and {9,pg) variables in ((37|l . and expressing one set in terms 
of the other, leads to 

a; = —2 In I — — - ) , p = —ipg tanff (39) 

\— mcosyy 

9 = arcsin {^—ie^^^p/rnj , pg = —me" cos (arcsin i^—ie~^^p/m)) 
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where some care is needed if cuts are encountered. These relations lead most easily to the result that p (t) lies 
on a complex circle, for real i;^ > 0, as a simple consequence of 9 (t) = 6 (0) + 2p0t, where Im 6 (t) = Im d (0) is 
an invariant for real E > 0. Comparing H39|l to H33|) . we see that a — exp {~2i9 (0)), so \a\ = exp (2Im0 (0)). 
Thus, X (t) turning points are encountered for real £^ > if and only if Im 9 = on the corresponding free 
particle trajectories. 

Actually, this canonical transformation works no matter what the energy is, even if it is complex, in 
which case the classical trajectories can have an elegant appearance. The Figure below shows a momentum 
trajectory for a; (0) = 1 and E = e^. The flow is between symmetrical complex fixed points at p± = ±e*/'^. 



The choice E — e^' would give a more symmetrical trajectory, with p{t) ~ —p{~t), and with a turning 
point. Of course, as indicated by the explicit transformation H39|) . complex energy would require extending 
the free particle trajectories into both the complex 9 and the complex pg planes. 

The classical trajectories and the classical canonical transformation which we have just described will 
appear explicitly in the quantum mechanics of the model, to be discussed in the next section, where they 
play pivotal roles in canonical integral representations of the wave functions and their duals. But before 
ending our discussion of classical dynamics and turning our attention to quantization, some comments on the 
canonical structure of this system would seem to be in order. Indeed, this issue is brought into particular 
focus by our previous remarks concerning classical canonical transformations. 

We might have proceeded as in the usual case of a particle on the plane, and assumed the real and 
imaginary components of a; = Re a; + ilmx and p = Hep + ilmp are independently canonical, up to some 
overall "isotropic" scale cr, with Poisson brackets (PBs) given by [Rea;,Rep]pg = cr, [Ima;,Imp]p3 — a, 
[Rea;,Imx]pg — 0, [Rea;,Imp]pg = 0, etc. Had we done so, we would then have had for any value of a 
the brackets [x,p]p-g = 0, [x,H]p^ = 0, [p, i?]pg = 0. However, this is not the canonical structure that 
leads to Hamiltonian dynamics as given in H^Of) . nor is it the structure that underlies the classical canonical 
transformation given in 1)37(1 . To obtain non-trivial equations of motion and other canonical transformations 
of variables, in this isotropic Poisson bracket approach, it would have been necessary to also conjugate some 
of the dynamical variables appearing in H and the generating function F. But we chose not to do this. 

On the other hand, if we were to suppose the brackets are different for the real and imaginary parts, and 
modify them to be anisotropic brackets (ABs) as 

[Rex, Rep] = cr + cos^ C , [Imx, ImpJ^g = cr — sin^ C (40) 
[Rea;,Imx]^g = 0, [Rea;,Imp]^g — 0, etc., then for any values of a and C wc would have 

[x,p]^^ = l, [x,H]^^^2p, [p, i7] - -2zm2e2- (41) 

This is exactly the dynamics specified in ((3()|l if wc identify df {x,p) /dt — [f (x,p) , H]j^^, in the usual 
way. If we were to go further and impose the additional condition that the bracket of x with the complex 
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conjugate of p vanish, then we would be led to 2a + cos 2C, = 0, hence a + cos^ C = 1/2 and a — so? C, = —1/2. 
This would give the most symmetrical form for the anisotropic brackets. For this choice the real and 
imaginary components would have brackets differing by a naive PT transformation, i.e. just a sign flip: 
[Rea::,Rep]^3 = 1/2, [Ima;,Imp]^g = —1/2. This is what we have done. 

We will briefly return to this issue at the end of the next section, in a quantum phase-space framework, 
by considering deformations of the anisotropic brackets. 



4 Quantum mechanics for V {x) = exp {2ix) 

Reconsider the PT symmetric Hamiltonian, H29|l . only now in a quantum setting. Again, granted that PT 
symmetric quantum theories make sense, this has to be the simplest periodic quantum potential problem 
one can imagine.^ The corresponding Schrodinger energy eigenvalue equation is 

'^+m'e^"^4'E^E^E (42) 

Since the potential is periodic, with period tt, Floquet's theorem assures us that we can always find a solution 
of the form 

^ {x) = e*^^0 {x) (43) 

where (a; + tt) = cf){x). That is to say, [x) has the form [x) — J^nei. CnC^™"^. These features are easily 
borne out, since the equation can be solved exactly. Change variables to 



Z = 7716*"" (44) 



Then -w^ — iz-§; and becomes 



az^ oz ' 

This of course is Bessel's equation, albeit in the perhaps unfamiliar situation where the variable z/m is on 
the unit circle. Series solutions arc explicitly given by 

±VE^ (-mV4)" 



n=o nil i 



n±^/E 



and exhibit features in accord with Floquet's theorem. For completeness, we also recall the elegant contour 
integral representation of Schlafli and Sonine ( 82 §6.2) which is valid for all ^/E and x. 

The contour begins at — c», with argw ~ — tt, proceeds below the real w axis towards the origin, loops in 
the positive, counterclockwise sense around the origin (hence the (0+) notation), and then continues above 
the real w axis back to —oo, with argw — +7r. 

The two choices of sign on the index of the Bessel function give independent solutions, so long as ^/E 
is not an integer. The prefactor in (|46() or H47() is just a free plane wave, for all real E > Q. The infinite 
sum (or Sonine's integral) provides periodic modulations and additional phases to the plane wave. Negative 
or complex E would give for the solutions either an exponential blow-up for x on the real line, or else 
non-periodic behavior for x G [0,27r7V], where in the latter case periodicity requires ^/E = k/N for some 
integer k. For these reasons negative or complex E is always ruled to be madmissible. In this regard, we 

^This example is mentioned in passing by Bender et al. without much discussion. It is considered in more detail 

in 1201 . Moreover, the model has been discussed extensively in the mathematical literature, at least as early as the 1950s 
EZIESlESlEnilSSllHniElEniEDEZI- The model with periodic boundary conditions on x G [0,27r] is an exception to the 
method used by Bender et al. to construct the dual space. We will discuss this fully in the following. 
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emphasize that our discussion appHes to x on the real hne, or else to 2; = me'^ lying on a closed contour 
surrounding the origin. For energy spectra corresponding to boundary conditions on other contours in the 
complex plane, see §5 in |2U| . 

So, the solutions of l|42|) for all real E >Q are (recall J„ (z) = (—1)" J-n [z) for integer index) 

, \ 1 / t / 7t\ t / itw f Jn (fne^^) for n even , 

/ N 1/t / ix\ T / ixw f for n even /.o1^ 

V',.(x)=2(Jyi,(me )-J_^(me ))^^=^ | J„ (me") forn odd ^^Sb) 

We have chosen linear combinations of J^.^ for later convenience. We have also chosen the phases of the 
eigenstates so that VT tpE i^) = ''Pe (^)- When E = 0, since "^b^o = 0' have a non-degenerate ground 
state, i^E^a (^) ~ ^0 (me") . This is unlike the case of real Liouvillc theory, for which there is no ground 
state 32 . 

In fact, note the lack of degeneracy whenever \/E is any integer. If \/E = n > 0, the two independent 



solutions of (|45l) are usually taken to be J ^{z) and lim Y^{z) = lim •^VEi^')'^°^{^^^^ ^ ,/e(^) ^ ^Yje 



latter being well-defined even as vi? ^ integer. But, although Yn yme^^^ is finite for real a;, nevertheless 
it can still be ruled out as an admissible solution because involves a logarithm multiplying a periodic 
function. Hence it has a periodically modulated term with an envelope linear in a;, so it is unbounded for 
— (X) < X <oo. Alternatively, Yn (me") is not a periodic function for x € [0, 27r]. So F„ (me") is not an 
admissible solution on the domains of interest to us. 

We plot (Re ^/^E (x) , Im'f/'E (2:)) and (ReV'l; (2^) ,Im?/;|; (a;)) parametrically, as functions of a;, for two 
cases. 




(Rct/ib (x) ,Im-0£;) and (Re?/'|;,Im'(/'|;), for E — Q (green and orange) and for E — 1/4 (blue and red) 
In the case £' = 0, it is sufficient to take x G [0,27r]. In the case E — 1/4, we have plotted ip^^^ and 
(^1/4) ^ ^ [— 27r, 27r]. The complete curve for the other wave function at this same energy, "01/4 ? is 
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obtained by rotating the complete "01/4 curve through 7r/2 about the origin, even though point-by-point, 
'01/4 i^) 7^ '3^*^^^V'i/4 i^)- Similarly the complete (^ipi/ij curve is obtained by rotating {^^1/4^ by tt. 



From these plots it is clear that (x) is not a bona fide probability density since it is complex. It is also 
more or 



less clear from the symmetry of V'^ that (0o)^ dx is real and positive, while /^jtt (^1/4) 
real and negative. These are the so-called PT inner products This leads to the problem of constructing 
the dual space for all the eigenfunctions. 

When \fE is not an integer, the "CPT" construction of the dual space 5, goes through. This follows 
because 



^%^{x)i>%^{x)dx^±b{E^~E2)y^2sxJT:^fE'^ , / ^+ (x) Jx) = (49) 
and sin ^tta/B^ on the RHS does not vanish when \/~E is not an integer. In this case, the dual functions 



can be identified with 

2 sin ( tt^/E 



and exhibit the expected sign flips for alternating bands of energies. However, when \ E \s an integer, 
sin {n^/E^ vanishes, the PT symmetric eigenfunctions are self-orthogonal, and the CPT construction of 

the dual space fails. This happens for the invariant subspace of non-degenerate eigenfunctions which are 
periodic on x G [0, 2it\. In the mathematical literature |42[ I8UI 177] these integer ^/E are known as spectral 
singularities^, a name seen to be appropriate from the behavior of the integrands for spectral resolutions of 
projectors. 

. . . ±1 



P± ^ , (x, y)= I 0± (x) (y) T^^dE (51) 

' ' Je^ 2sinf7r%/;Ej 

Henceforth we consider only the subspace of 2tt -periodic eigenfunctions. 

A remarkable 27r-periodic biorthogonal system Elements of the dual space for the 27r-periodic eigen- 
functions are given by the Neumann polynomials as presented earlier. The key orthogonality relation is 
now 

1 f^"" 

— Ak (me") J„ (me") dx = 6kn (52) 
Jo 

where A-n (z) = (—1)" An (z). This gives a remarkable biorthogonal system whose properties warrant more 
discussion. The CPT method of constructing the dual space, i.e. simply changing normalizations and phases 
of the PT transformed functions, does not provide a bi-orthonormalizable set of functions in this case, since 



1 



Jfe (me") Jn (me") dx = Sko6„o (53) 



This follows because the Js are series in only positive powers of e". So all the 27r-periodic energy eigen- 
functions are self-orthogonal except for the ground state. 

The effects of the Hamiltonian on the dual space The An obey an inhomogeneous equation |82j . 

2nme" for n odd , . 

2m^e^" for n even ^ ^ 



-^An (me") + (m^e^" - n^) An (me") - | 



This leads to an interesting variant on the usual method of showing that eigenfunctions and their duals are 
orthogonal for different energies. Consider J^^ Ak (me") HJn (me") dx = J^^ Ak (me") J„ (me") dx. 

^For other examples of this effect in PT symmetric theories, see 1741 . 
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Integrating by parts gives 

r.2ir 



/ Ak (me") HJn (me") dx = J, 
JO Jo 



me") i/Afc (me") 

2Tr 



^27r C 2fcm e"J„ (me") da; for fc odd 

= /cM J„ (me") Afc (me") dx + <^ (55) 
" [ 2m^ Jp'^e^"J„ (me") dx for k even 

There are no boundary contributions upon integrating by parts since all functions are 27r-periodic. More- 
over, the inhomogeneous terms are orthogonal to the space of energy eigenfunctions. That is to say, 
J^^ e^^Jn (me") dx = = J^^ e^^^ Jn (me") dx since the Js are series in only positive powers of e". So 
even though (|54() is inhomogeneous, we still arrive at the usual conclusion. 

(n^ - fc^) / Ak (me") J„ (me") dx = (56) 
Jo 

hence A^ (me") J„ (me") dx = if /c^ 7^ rt^. Recall that (z) = (—1)" An {z). 

Integral representations for E = and quantum equivalence to a free particle on a circle The 

27r-periodic Bessel functions are in fact the canonical integral transforms of free plane waves on a circle, as 
constructed in this special situation just by exponentiating the classical generating function presented earlier 
in 1(23. Exphcitly, 

J„(me")==— / e^'^^'^'^^e-^^dei , nel (57) 

^TJ" Jo 

The essence of the quantum equivalence lies in the identical action of H and i/froo — —dg on the periodic 
kernel. 

^^ime" sine ^ ^^^^^g^me^ sin S ^ ^^2^2ix ^ [m^^^ ^ gime" sin 9 (-gg^ 

While (|57(l is a well-known integral representation of the Bessel function "F? , the interpretation as a canonical 
transformation is usually overlooked in most mathematical physics texts. That is unfortunate. This is 
an extremely straightforward and transparent connection between the quantum and classical theories that 
serves to distinguish this model. (It doesn't get any better than this!"^) 

The integral transform is a two-to-one map from the space of free particle plane waves to Bessel functions 
(e^™^ — > (±1)" J„ upon using J_„ = (— )" J„), just as happens in the real Liouville case {^e^"^^^ Kik) 
[13 US], although here we have selected normalizations for the wave functions so that the transformation 
uses the same normalization magnitude for all n (with at most only an n-dependent phase), unlike the real 
Liouville integral tranformation, where the magnitude of the normalization varies with k in order to obtain a 
conventionally normalized continuum wave function from the MacDonald function Kik (e.g. |26| . 5th eqn.). 

When restricted to the combinations e™* -I- (— )"e~™^ the map is one-to-one and invertible. On the 
other hand, acting on e*"^ — (— )"e~*"^ the map projects to zero. When restricted either to the subspace 
of right-moving plane waves, or to the subspace of left-moving plane waves, the map is also one-to-one, but 
not invertible to the same subspace of plane waves. Inverting this transformation on J„ by using the same 
kernel always gives a combination of left- and right-movers. (See the discussion and Table below.) This 
left-right combination is unavoidable, given the properties of the kernel exp {ime" sin^) and the symmetry 
J_„ = (— )" Jn, although the combination of plane waves can be changed by choosing a different kernel in 
the integral representation of J„. (See the discussion following the Table.). In any case, combinations of 
left- and right-movers result from inverting the map using the kernel we have chosen. 

We note the effects of index raising/lowering operations on the complete integrand in 



^The same situation applies to the real Liouville model. See [TTI ITSl HSl 1^ l?7l 1^ . 



12 



Thus, acting on the integral representation (|57|l and integrating away the 9 derivative term, we obtain the 
usual recursion relations only expressed for Bessel functions on the circle. 

T Jn (™e") = Tme*" J„±i (me*-) (60) 

We also note in passing that the exponentiated classical generating function has a well-known interpretation 
as an element of the Poincare group on the plane ( 81 Chapter IV). 



Inverting the transformations and the role of other classical trajectories in the quantized model 

The canonical integral transformation from Neumann polynomials to plane waves using the same kernel as 
in H57(l warrants discussion. It is simply 

This follows from the biorthogonality of {Jj,Ak], (|52|l . and the well-known generating function (see 
9.1.41) 

oo oo 

gi.sine^ ^ J„(z)e^"^ = Jo(z) + ^ J„(z)(e^"^ + (-l)"e-™«) (62) 

n — — C30 n— 1 

as well as J-_„ (z) = (— 1)" Jn(z). The transform H52|l involves integration around the circle, z = me^^ , 
X e [0, 27r], of a Neumann polynomial times the canonical transformation kernel. But, were we to make 
slight deformations of the circular contour (as allowed by the analytic form of the Neumann polynomials and 
the kernel in the neighborhood of the circle) the result of the transformation would be the same. Therefore, 
the transformation may be carried out just by integrating along a closed contour in the z — me" plane given 
exactly by following a single classical Liouville particle x-trajectory, for E > Q. 

To invert the transformation and go from plane waves to Neumann polynomials using the same kernel 
ultimately requires taking free particle plane waves off the circle and integrating up the imaginary axis. We 
may do this by first integrating along an arc of the circle to the point 9 = —a, and then by following the path 
of a classical free particle trajectory for Ivcipg > 0, Repe = which evolves to ioo — a a,s t —^ oo. Explicitly 
(IH2 §9.14 Eqn(3)) 

1 noo-\-ia 

An (z) = -Snz / (e"^ + (-l)"e-"'') e""'''"'^'' coshMi^ (63) 

2 Jo 

where jo; -I- argz| < 7r/2. Integrating by parts picks up a contribution at when n is even, but not when n 
is odd, resulting in Schlafli's version of the Neumann polynomials (see §9.3-34 in |H21)- Exphcitly 

n / 1 \ i-oo+ia 

Sn {z) = -( -An {z) - cos2 (n7r/2) = / (e"" - (-)" e""") e-^-^^^d,? (64) 
n \en J Jo 

Alternatively, in terms of the original free particle variable 9 = id. 

9/1 \ fioo-a 

Sn iz)^-{ —An {z) - cos2 (n^/2) = i (-)" / (e™^ - (-)" e-*"«) e*^-"«d0 



n \s. 



zoo — a 



An (z) = En cos' (n7r/2) + zsn (-1)" - / {e'"' - (-)" e-™^) . 

^ Jo 

(65) 

cos (n9) d9 for n odd 

At least three wedges (choices for a) are required for an open cover of a full circle for z, although two wedges 
will cover all but two points on the z-circle. 

We give a summary Table of the canonical integral transforms using the kernel exp {ime^^ sin 0) . The 
paths of integration are described in terms of relevant classical trajectories. In a well-defined sense as 
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indicated in the Table, only integrations over particular, representative, classical trajectories are required to 
change wave functions for the free particle into those for the Liouville system, and vice versa. 



Integral Representations 



with 



Relevant Trajectories 



J„ (z) - — / e*"^'"« (e-™^ + (-1)" e™'') dO 



z = me 



on Ke pe > , Impg = 
free particle trajectory 



icosfl / e 

' — oo 



V„ (z) = e"'« + (-)" e 



9 real, 
sin 6* ^ ±1, 



X on ReE = = ImE 
X (0) and tt 
Liouville trajectories 



(■2) - ^ / e^'""^ (e"™^ + (-1)" e™«) cos0d0 
2* Jo 



la + argz| < 7r/2 
i.e. |a + a;| < 7r/2 



61 on Repe < , Impe = , 6i (0) = 
andRepe = , Impg > , (0) = a 
free particle trajectories 



— /e*^^'"''A„ (z) — = (e™« + (-)" e-™«) ^ 
27ri / z ^ '2 



a; on Re£; > , ImE = 
Liouville trajectory 



Analyticity is in control here, in the Table. Up to some clearly discernible limitations, variations of the 
integration paths will not change the results. Indeed, some deformation of the integration path may be 
necessary to have it exactly coincide with the indicated classical trajectory (for example, in the last case). 

As previously remarked, the first and last rows of the Table are immediate consequences of the biorthog- 
onality of {Jj,Ak}, (|52l) . the well-known generating function, H()2|) . and J_„ (z) = (—1)" Jn {z) as well as 
A-n {z) = (— 1)" A„ (z). The second row follows from the first row and the Fourier integral representa- 
tion of the Dirac delta. This particular relation is sometimes written an as integral representation for the 
Chebyshev polynomials, after translating 6 9 + 71/2 and defining u = cos 9, whereupon it becomes (see ^ 
11.4.24) 

cos {n9) EE r„ {lu) = -i'Vl -t^^ / e"'^" J„ (s) ds for < ^ (-gg) 

The third row of the Table is obtained just by substituting 9 — ii!} in lt)3|l . with no further insight on our 
part, but it may well have a more physical interpretation |46| . 

For the kernel being used, as presented in the Table, the map is {J„} ^ {e"*^ + (— )" e~'"^} ^ {^n}- 
The intermediate subspace of plane waves is not chiral. However, projections onto chiral subspaces may 
also be constructed. 

Chiral kernels The result (|52|l is a bilinear but non-symmetric kernel that maps linear combinations of 
plane waves onto integer Bessel functions. Alternatively, the kernel maps Neumann polynomials onto equal 
mixtures of left- and right-moving plane waves. But obviously, we may project onto half the space of plane 
waves in an infinite number of other ways, say by choosing to keep in the sum over n only one of exp (±in9), 
for all n > 0, and thereby produce other maps that act nontrivially on the complementary set of plane 
waves. 

For example, the left-moving chiral subspace of plane waves is mapped one-to-one onto Bessel functions 
by the chiral kernel'* 

C30 

5(z,0) = ^ J„(z)e™'' (67) 

n=0 

Thus S : {e~™^ | n > 0} ^ { J„ (z)} . Alternatively, all Neumann polynomials are "right-mapped" by the 

''The right-moving chiral subspace can be similarly mapped onto Bessels using S (z, 9) = S {z, —0). Also note the radius of 
the free particle circle can be modified, e'^ — » Me*^ as in 2 = me'^, by giving 9 a constant imaginary part. 
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kernel into the space of right-moving plane waves, {An (z)} : S — > {e™^ | n > O}. Explicitly 

1 /-Ztt 1 I- J 

J„(^)=2^y^ S{z,0)e'^"'de , e"'' ^—j-A^(z)S{z,0) (68) 

Moreover, note the null projective property on almost all of the right-moving subspace, S : {e™* | n > O} ^ 
{0}, as well as on all Bessel functions except for Jq, {J„ (z) \ n > 0} : S ^ {0}. 

The chiral kernel can be expressed as special cases of Lommel's functions of two variables (|H21 §16.5-16.59) 

S {z, 0) = C/o {ie^'^z, z) - iUi {ie^'^z, z) (69) 

but regrettably, there does not seem to be an elementary closed-form for 5*. From rewriting (|62|l as 
S^-oo ^""^n (^) = ~ 7)1 however, an integral representation follows immediately. 

S{z,e) = ^l e^<*'i)-l-dt (70) 



27rij|t|>i 

Here the t contour is counterclockwise around the origin with \t\ > 1.^ The simple appearance of H7U|I is 
somewhat deceptive. 

As in the case of the previous kernel, e'^ ^ , this S establishes a quantum equivalence between complex 
Liouville and free particle systems. Only now the equivalence is with free chiral particles on a circle. Once 
again, this is directly evident from the identical action of H and -fffroc = Pg on the kernel. 

00 

HS (z, 6) = ffft.ec5 (z, e)^Yl (■^) (71) 

n=0 

This may be checked against the partial differential equations satisfied by the Lommel functions in l|69|) . or 
the contour representation in (|70|l . From either of those results, or by manipulating the sum in (|7T|) . there 
follow other forms for this bilocal Hamiltonian kernel acting on the two spaces. 

The chiral kernel can be inverted on the space of Bessel functions by the formal kernel S^^ (0, z) . 

(0, z) = Y, e-"''An (z) , e-"^ ^ i^S-\0,z) J„ (z) (72) 

n=0 ^ 

Thus S-^ : {Jn (z)} {e-™^ | " > O}, and therefore S'^S = 1 : {e-"'^ | n > O} {e"'"^ | n > O} and 
SS-^ = 1 : { J„ [z)} {Jn (z)}. Alternatively, {e™* \ n > 0} : {An (z)} since formally 

An (z) = r e"''S-^ {0, z) d0 (73) 

The Liouville and free Hamiltonians act differently on S^^ {&iz) due to the inhomogeneous terms in H54|l . 
Acting term-by-term and summing the formal series we find 



{H - i7ft.ee) S-^ {0, z) = 2zide ^ e""^ -f 2z'^ ^ 

2i0\ 

r2 (74) 



e 

odd n >0 even n >0 

2ze-*'' (1 + e-2*«) + 2z2 (l - e^^''^) 



(l-e-2»») 

This has singularities as a function of 0, but these can be avoided in any integrations by changing the radius 
of the free particle circle (i.e. by taking Im6' ^ 0.) The RHS of (|74|l is also manifestly orthogonal to the 
right sector upon integration 

This formal inverse allows us to express the equivalence to a free chiral particle in a somewhat heuristic, 
but very concise way. 

-fffrcc chiral = S HS (75) 
^If \t\ < 1 for the contour, instead of S it yields I — S {z,8 + it). 
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Or more explicitly 



fffree chiral (0', &) = 7^. f S"' {9', z) HS {z, 0) = V e""^ n^e™^ = (76) 



The last expression here is the usual form for the free chiral Hamiltonian kernel, and it is manifestly not 
local, whereas the complex H is local. 

As a technical point, we note that {9, z) is indeed a formal kernel, since it is a divergent series, given 
the normalizations of the Neumann polynomials. Nevertheless, it can be Borel summed in closed form (a 
result due to Kapteyn, see §9.16). 



oo ,.00 



n=0 







where the integral is recognizable as the exponential integral function Ei. Conversely, the sum in (|77|l 
is the asymptotic expansion of the integral, the latter being well-defined and obviously convergent when 
Re (1 — t^) z/t < 0. As generally remarked following (|i;-{|l . the convergence of this generating function for 
the dual polynomials can be improved by making compensating changes in the normalizations: {Jj, Ak} —^ 
{ZjJj, Ak/Zk}- Thus Zn can be chosen to make X]^o F"^" ('^) /■^" absolutely, uniformly, and rapidly 
convergent. The price to be paid is that J2'^=o ^nJn {z) e™^ will then be a divergent series.^ 

Other non-symmetric kernels may also be constructed. For simple examples, reconsider those in (123(1 and 
((25|l . Of greater interest, perhaps, are symmetric kernels given as sums of bilinears where both functions 
are of the same type. It is often possible to express these in closed form as well. 

A simple kernel for the bilocal norm While it may be difficult to find a closed form for the kernel 
K (x, y), since our normalizations for the An (me*^) make this a divergent series, it is a simple task to find a 
closed form for J {x, y), the inverse of K on the spaces of interest. We recall Neumann's addition theorem, 
l|27|l and (|28f) . only this time we choose more convenient normalizations to write 

00 

J (x, y) = Jo (me-'" - me'y) = ^ e„ J„ (me"*") J„ (me*^) (78) 

n=0 

where Eq = I , Sn>i — 2. Formally, the determinant of this kernel is 

det (J) = exp (Trln J) = exp / ln{J{x,x))dx 

Jo 

= exp / In (Jo (2jmsin2;)) dx 
Jo 

p27T 

= exp / In (/q (2TOsinx)) dx (79) 
Jo 

where Jo is the modified Bessel function (see (|93f) below). Thus det ( J) never vanishes, for any real m, 
implying the existence of the formal inverse for real m and x € [0, 27r]. With the normalizations in (|78|l that 

''The most "balanced" asymptotic index behavior for rescaled eigenfunctions and their duals is achieved as follows. 

n!2" J„ (z) z" , ^— A„ (t) — 

Then the leading behavior of the terms in the bilinear sums, modified accordingly, is 

-.00 -. 00 -. 

n = JV ^ ' ^ ' n=JV st 
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formal inverse is 



K (x, y) = Y^ —An (me") An (me"'^) (80) 



n=0 



N 



These exact results should be compared to those for the spectral resolvent, in models where V {x) = (ix) 

As in H13|l . the manifestly hermitian kernel J {x,y) can be used to evaluate the norms of arbitrary 27r- 
periodic states 

oo 

(a;) = ^ Cny/e^Jn (me") (81) 
through use of the corresponding dual functions 

oo 

V'dual {x) = ^ c*,An (me") /Ve^ (82) 

n=0 



— 2 / dx dy ^duai (a;) J {x,y) V^duai (y) = ^ |c„|^ (83) 

"'0 "'0 .,_n 



IIV-II" 

(27r) .u „^o 

We note that J{x,y) converts the dual polynomials into conjugated Bessel functions, and conjugated poly- 
nomials into BesselsJ 

I /■27r 



J: (me") = J„ (me-") = / dy J {x, y) A^ (me'«) 

^Tl'Sn Jo 



^ /.2ir 



J„ (me^«) = / da: A„ (me"") J (x, (84) 



27re„ Jo 

Hamiltonian kernel on the dual space We also note that J (a:, y) converts H into H* to produce a 
bilocal Hamiltonian kernel equivalent to either, in the following sense. 

-Q^+ "I'e-^" j J {x, y)^i^-—+ m^e^^yj J {x, y) = E ^'^«"^» {me-'^) J„ (me*^) (85) 



The sum converges absolutely, uniformly, and very rapidly. Evaluating this sum in closed form produces an 
explicit hermitian kernel. 

H (x, y) = H* {y, x) = m^l— -—- = ^ n'enJn (me"") J„ (me'«) (86) 

n=0 

As a series in m (e^" — e*^) this bilocal Hamiltonian is rapidly convergent 

H ;!!e---.»f: r-f)'" (87) 

^ '^^ 2 ^ 4"n! n + 1 ! ^ ' 

n— ^ ' 

and as a series in e"" or e*^, it only involves non- negative powers of either. This Hamiltonian kernel acts 
on the dual space to give energy operator matrix elements, or averages as in (|15|) . 

(H)^ = ,2,1 ,,|2 / dx / dy V^duai (a:) H {x, y) ^duai (y) = E / H^H' (^8) 

(27r) IIV-II Jo Jo „^o 



It is sometimes convenient to re-express this in terms of the action of H (x, y) on the dual polynomials, 
similar to 



(me") = nV„ (me'") = / dy H (x, y) ^„ (me^^) (89) 



'^Simply changing the sign of x in J{x,y) yields a non-hermitian kernel that converts A„ (^rne^y) into non-conjugated 
Jn (me'^). 
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Another map between H and H* on 27r-periodic functions Consider [201 integral kernels depending 
on a parameter s. In particular, consider the very simple kernel 

K {x,y;s) = eyi-pk{x,y]s) (90) 

k [x, 2/; s) = ^ Qe^(-"+^) - se*(-=^-^) - se'(-+^)^ 

k (x, ?/; s) = k* {y, x; s) for s e R 

All terms in these expressions are periodic in x and y, separately, with period 2tt. Once again we have 



dx^ 



+ m^e-^"^ K [x, y; s) = (^-^ + m^e^^^^ K (x, y; s) 



(91) 



So this kernel also accomplishes the task of converting H into H* and vice versa. Using this kernel, we find 
the integral transforms* 

J: (me*^) = J„ (me-*^) = — f-^ / if (x, y; s) J„ (me^«) dy 

J„ (me'^) = — / J„ {me-'-) K (x, y; s) dx (92) 

' 27r /„ (ms) Jo 

which should be compared to (|84ll . For real ms 7^ 0, the modified Bessel /„ (ms) is real and it has no zeroes, 
as is evident in 

^''("^■^)-(-) g fe!r(i + , + fc) (-) 

For all real m and s, if (x, y; s) is an hermitian kernel, K (x, y; s) = K [y, x; s)*, and it provides an explicit, 
one-to-one, invertible, closed-form transformation between the cigcnfunctions of H and H* . 

Another way to state the result H92I) is as an expansion of the kernel in terms of Bessel/Neumann 
trilinears. (A similar relation holds for real Liouville theory. See 1261.1 From the transforms given, and 
the completeness of the J„ (z) on positive powers of z, we deduce that 

K (x, y; s) = Jq (me^") Iq (ms) Jq (me*^) 

00 

+ J2 (-1)" Jn (me-") /„ (ms) A„ (me^^) + (-1)" A„ (me"") /„ (ms) J„ (me^^) 



00 



+ J2 J J {me-'-) cjk is) Jk {me'y) (94) 

where we have not determined the coefficients Cjk (s) (except for the omitted cqo = = cqi = cio)- 

More information may be incorporated into the sums expressing K (x, y; s) by imposing total symmetry 
among e"", e'**, and is = e'^. Since J„ (me*^ = ims) = i"/„ (ms = —ime'^) we have 

J„ (me") J„ (me*^) = / ^ (2:, y, z) J„ (me^^) dz (95a) 

27ri" /n 



S" (x,y,z) = exp 



+V-^) _|_ gi(a;-y+2) _|_ ^i{-x+y+z) 



K {~x,y;s = -ie") (95b) 



where 5 is unchanged by permutations of its arguments. It follows as above that this symmetric kernel is 

*A similar kernel may be used in real Liouville theory to effect an integral transformation from one K,y to another, leading 
to MacDonald's identity (l^'M §13.71). See I2til and references therein. We were not able to find I92i in the literature. But it 
is not too difficult to work it out. 
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of the form 

oo 

S {x, y, z) = Jo (me") Jq (me*^) Jo (™e*') + ^ c,fe„ J, (me") Jk (me*^) J„ (me*^) (96) 

j.k,n—l 

oo 

+ ^ i"A„ (me") J„ (me^^) J„ (me*^) + J„ (me") A„ (me^^) J„ (me*^) + J„ (me") J„ (me*^) A„ (me*^) 

n=l 

We have not determined the symmetric coefficients Cjkn- Regardless of their values, they drop out of (|95a|l . 
Their relation to the previous coefficients is 



Cjk (s) = Sjki''Ak (ims) + ^ cjkni^'In {ms) (97) 

n=l 

The Cjfen could be determined by evaluating J^^ S (x, y, z) An (me*^) dz through a direct calculation. 

Classical behavior of time-dependent expectation values As a final comment on the structure of 
this model, we remark that it is straightforward to take linear combinations of energy eigenstates to produce 
near-classical time-dependent behavior for expectation values, including those which may be complex, such 
as 

(P)^ (0 = . 2 / dy V'duai {y, t) i-idyiP {y- 1)) (98) 

27r /.2ir 



27r||V||^ Jo 



dx dy Tl>duiii{x;t) {-idyJ {x,y))il!duiii{y;t) (99) 
Jo 

Here we have defined the time-dependent function and its dual as expected from (|81ll and H82|) . 

{x; t) = Yl Cne""'* V^Jn (me") , ^duai {x; t) = Yl <e™'*A„ (me") (100) 



n=0 n=0 



We have also used 'tp{y;t) — Jq^ dx ^duai {x; t) J {x , y) to write the expectation value in terms of the 
non-hermitian bifocal kernel —idyJ {x,y). 

In fact, a simple combination of only two energy eigenstates will give (p)^ (t) that evolves around a circle 
in the complex plane, whose center is on the real axis, with constant angular velocity around the circle. We 
believe that only slightly more complicated coherent states [221 can be constructed to reproduce the exact 
time-dependence evident in the classical motions, where the angular velocity around the momentum circle 
is not constant. All this is in accord with the time derivative 

-Ti ip)i, (t) = -i (b, H])^ [t) = I dx V-duai (a;; t) [p, H] {x; t) 

"■'^ 2t:\\iP\\ Jo 

^-2^m^e'^^)^it) (101) 

Classical evolution of (p) ^ (t) follows if the state V' is such that the last RHS evolves classically. Since we 
also have 

I (e^")^ (t) = -^{F-,H])^ (t) = 2^{{p,e'-})^ (t) (102) 

classical behavior is guaranteed to the extent that ({p, e^"})^ (t) coincides with 2 {p)^ (t) (e^")^ (t). (This 
is just Ehrenfest's theorem for this biorthogonal system.) In this regard, it is crucial that we are able to 
replace by HAn under the integrations over x in H101|l and l|102(l . Fortunately, the inhomogeneous 

terms in (|54|l do not contribute to the expectations under consideration, so the replacement is permitted 
among the integrands in these cases. 

Remnants of classical behavior or the classical limits of quantum systems are often most clearly under- 
stood from the point of view of phase-space quantum mechanics, i.e. deformation quantization [86) . In that 
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framework, let us briefly return to the issue of the anisotropic brackets. The classical anisotropic brackets 
described earlier are just h limits of Moyal brackets built from anisotropic Groenewold products acting 
on phase-space distributions, where the phase-space variables are allowed to be complex. Thus 

[f {x,p) ,g{x,p)]j^j^ = lim (103) 
n^O III 

where the associative but non-commutative anisotropic product is given in its most symmetrical form by 

^ d Be X d Jlc p— d U.C p d U.C d Im j: d Im p— d Im p d Im (^104) 

Clearly then, = xp + ih/2, p'^x = xp — ih/2, and = 1 as in (|41|l . as well as [Rex, RepJ^g = 1/2, 

[Imx,Imp]^g — —1/2, etc. Therefore the previous continuation of the classical dynamics is recovered in 
the h —^ limit. In fact, the quantization of the system can be carried out in full detail on the complex 
phase-space through the use of this star product acting on Wigner functions built from bilinears in the wave 
functions and their duals. (A thorough discussion of this approach to quantization is in preparation |28|.) 



5 Magnetic field effects 

Magnetic fields can produce some remarkable changes in the structure of a complex biorthogonal system, 
even in the simple physical set-up where charged particles move around a solenoid containing magnetic flux 
$ with only a constant vector potential A cx $ on the surface of the solenoid, in addition to the complex 
potential energy imagined to be already present.^ As an example, we add such magnetic effects minimally 
to the previous Hamiltonian, to obtain the form (v = —eA/c) 

H ^ {p + lyf + m\^" (105) 

Here p is the momentum canonically conjugate to x, so classically there is no discernible change from the 
case 1/ = 0. Classically, we have dx/dt — [x, H] = 2 {p + v), so H = ^ [dx/dt)^ + m^e^", and for a given 
energy the trajectories are the same as obtained in the u = Q case. A corresponding classical canonical 
transformation taking H Hi^cc = Pg is now given by 

d d 
F = -i^x + me" sind , p = —F = -v + ime'''" s\n6 , pg = -—F = -me" cos6 (106) 

ox do 

Thus there are no significant changes, classically. 

But consider the QM of this model, and consider only 27r-periodic solutions, since the effects of 7^ 
are most dramatic in this situation. As an operator 

H = {-id^ + uf + m^e^" (107) 

and for generic v the set of periodic eigenfunctions now consists of two sectors given by 

U e"*"^" J,±„ (me") = {v ± nf g-*"^" J,±„ (me") (108) 

for n G Z>o. These eigenfunctions are all even under PT, for real v. For want of better terms, we will call 
these the right and left sectors since for generic v all right-moving and left-moving plane waves are contained 
in the spans of {e~^'^'^ J^+n \ n S Z>o} and {^e~'''^^ J^-n \ n G Z>o}, respectively. ^° When v ^l, Jv±n give 
two independent states for all n > 0. Thus the number of states has doubled from the 1^ — situation, 
no matter how small 7^ is.^^ The effect of ^ Z is to restore all left-moving plane waves to the span 
of the energy eigenfunctions, in sharp contrast to ly ~ where only the right-moving plane waves could be 
expanded in J„ (me"). 

^Alternatively, non-hermitian H could be obtained just by considering complex vector potentials for a particle on the plane, 
without any exponential interactions. See 1391 . 

^"Yes, n = is in both sectors, as defined. If this is troubling, then feel free to define three sectors to avoid any overlap. 

^'^This 2 — > 1 accounting of the states, as ^ 0, is complementary to the 1 — » 2 state counting for linear "wall" potentials, 
as the slope of the walls goes to zero. 
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So the energies are given on the two sectors by 

E±n = i'^±nf for n = 0, 1,--- (109) 

For generic v, E±n are not degenerate. But if — 1/2 = A e Z, eigenvalues are doubly degenerate. For 
example, if i/ = 1/2, the degenerate pairs of energies are -Eq.-Ii -Ei__2, etc. In this case there are two 
independent states for each energy. 

Otherwise, when i/ G Z, most if not all of the J^-n states meld into the J^^n states^^, since the 
Bessels with negative integer index are not independent of those with positive integer index, i.e. Ju-n {z) = 
(— 1)"^"" J„_i^ (z), for n and v G Z. Sorting out the independent functions leaves just {J\ \ X G l^^-y]. 
For instance, when v G Z>o the independent states are e"*'^^ Ji^_„ for — n S {0, 1, • • • , i^} (alternatively 
n G {i^jiy — 1, ■ ■ • ,1, 0}) and e"'"^ J,y+„ for n G Z>i. That is to say, when v G Z>o, the independent states 
are nondegenerate with energies 0, 1, 4, 9, • • • , exactly as in the v — situation (although the wave functions 
are different for v = and v 0). Moreover, when G Z>o only a finite number of left-moving plane 
waves are in the span of the energy eigenfunctions, namely: {e~"^^, e"*^""^-*^, • • • ,e~"}. Similar results 
apply when v G Z<o, although in this case only right-moving plane waves are spanned by the eigenfunctions, 
namely: ^gA'^l^ ^ g'-i'^+Wl)^ ^ Q^i^+H)^ ■ ■ ■ }; with the first right-movers absent from the span. But again, 
the energies are 0, 1, 4, 9, • • • , exactly as in the v — situation. In any case, all of these situations involving 
integer v are (gauge) equivalent through energy preserving one-to-one maps connecting the various Bessel 
functions. Map details are left to the reader to establish. 

For generic v, the dual space may be constructed for all eigenfunctions, but again not by the naive 
FT method. Every function in the right sector, e~*'^^Ji,+„ (me'^), is once again self-orthogonal (except 
n = 0) as is every function in the left sector. So, an appropriate dual function is not proportional to the 
eigenfunction, rather it is given by e^^^J-^-n {me") /C„ where Cn — ^ J^^ J-u-n {me") Ju+n {me") dx. 
This dual function is not an eigenfunction of H . Instead, it is an eigenfunction of 

H ^ {-id^ - vf + m^e^" (110) 

obtained by flipping the sign of the vector potential^^, with the same eigenvalue as the associated right-sector 
function, 

H e'"^" J_,_„ {me") ^ + nf e*"^" J_,_„ {me") (111) 

Similarly the dual for every left-sector eigenfunction, say e~*'^^Ji,_„ {me"), is e*''^J_iy+„ {me") /Z?„ where 

Dn = ^ J-.+n {me") J,_„ {me") dx, and H e'''" J_,+„ {me") ^ {v ~ nf e"'''J-,+n {me"). The 
duals for right- and left-sector eigenfunctions are guaranteed to be orthogonal to left- and right-sector 
eigenfunctions, respectively, for generic v where E±n are not degenerate, by standard arguments (as in (|54|l 
et seq.). But ii f — 1/2 = X G Z, further consideration is needed, since right- and left-sector eigenvalues are 
degenerate. An adaptation to biorthogonal systems of the usual degenerate state methods in conventional 
QM suffices to handle this situation for both right and left sectors, in a straightforward way. 

Canonical integral representations for the wave functions and their duals may also be exhibited for this 
model, as given for example by (|47|l with ^/E replaced by ± n and —v ± n. 

Rather than dwell on this, for simplicity we now limit our discussion to the right sector for any // ^ Z<o. 
This is an invariant subspace spanned by {e~^'^^ J,^+n {me") \ n = 0, 1, • • • }, and as previously pointed out, 
all right-moving plane waves are contained in this subspace for v ^ Z<o. As an alternative to calling 
it the right sector ,^it is the analytic sector in the variable z — me". For this subspace, rather than as 
eigenfunctions of H, the duals may be constructed as finite polynomials in 1/z in a manner that closely 
parallels our previous analysis of the 1^ = situation for 27r-periodic functions. 

Gegenbauer's generalization of the Neumann polynomials and the corresponding construction 
of the dual space In this generalization, the index n is always > and integral, while the index v is 

^^That is to say, integer u are exceptional spectral points. See I5UI 1471 . 

^^Note this is not what would happen under CPT, where "C" here really is charge conjugation, since u oc eA and this does 
not change sign under CPT. 
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allowed any value except —1, —2, —3, • • • . Negative v are allowed, just so long as they are not integral. An 
explicit definition is §9.2) 

^ 2''+"(i/ + n) r(i/ + n-fc) /zi;\2fc 

H - ^1.1 E ^ ( 2 ) (112) 

fe=0 



-u;"+i V {v + 7i-l) 2\{v + n-l){v + n-2) 

The Neumann polynomials are a special case: An^^^o{w) ~ £„0„(ii;). (NB A„^,y=o ("u^) — '^■^n{w) in 
owr previous notation. Sorry about that!) The Fourier transform on a simple counter-clockwise contour 
enclosing the origin once is 

1 °° 

—§e^^^A^^Aw)dw^2^i^{v + n)T{v)C'^{z) where ^ ^ rc„^ (z) (113) 

^ n— 

That is to say, the contour Fourier transforms of Gegenbauer's polynomials An,i, are just the Gegenbauer 
polynomials C^. 

The elementary Cauchy kernel is now given by 

— = ^ A„,, H z- (z) for \z\<\w\ (114) 



w — z 



This follows from expanding in powers of z and writing z^ for non-negative-integer fj, in terms of Bessels 
(|82j §5.2). 

„ A (ii + 2n) 2^r (u + n) , , , , , 

^^^y-U^ ^J^+2„ z for M^Z<o 115 

n=0 

Other kernels may also be constructed in closed form (for example, see H123II and H124(l below). 

These polynomials may be used to construct the dual space for the right sector. To this end we define 
for all n G Z>o 

^—iiyx 

V'n,^ (me") = Zn—^J^^+n (me") 

Xn,u (me") = -—A„^, (me") (116) 
where Zn is a normalization to be chosen for convenience, later (see H121f) '). As explicit series 

k=0 ^ ' 

, 2'^+"(i. + n) _„„^^^ r(i.-|-n-fc) /m\2fc 

A:=0 

Both sets of functions are clearly 27r-periodic. The energy eigenfunctions obey Schrodinger's equation 

HxPr^,. (me") = (^-i^ + 7/.„,, (me") -f m^e^"^™,. (me") = (i^ + nf i/;„,, (me") (118) 
while the finite polynomials obey the inhomogcneous differential equation 

(^H-{v + n)2) x„,. (me") = (^-i-^ - Xn.. (me") + {m^e^" - [y + nf^Xn,. (me") 

n odd 

(119) 



f £(!±|(!±1)) X 2me" for 

_ 2-{v + n) 

r(f+|n) 2 2?"t r 

—7 — i—i- X m^e for n even 
r(i+^n) 
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The two sets of functions satisfy the biorthonormaUty condition 



1 

2^ 



XkM (me") ipn^^ (me") dx = 4n for fc, n G Z>o (120) 



The proof for fc = n is obvious from the series (|117|l . The proof for k ^ n proceeds by first acting with 
H on V'njiy under the integral, then integrating by parts to act with H on Xk,v^ arid finally comparing the 
results, i.e. just following the same line of argument as in the case of a single exponential potential, H54|l et 
seq. But note the importance of the v sign change for H compared to H . 

Thus for any v ^ Z<o we have obtained for the right sector a biorthogonal system that is manifestly 
in one-to-one correspondence to that of the v — situation: {xk.vT'^n.v} < — * {^fc,>/n}- We leave the 
construction of explicit integral transformations between the v ^ Q and the = systems as an exercise. 
(Hint: see 12 §11.3) 

On the dual space of the right sector, the bilocal norm kernel can once again be obtained in closed form 
if we choose 



Then the hermitian kernel is 



^2 _ V4^ {u + n)T {2v + n) 
" " 2-r (i. + 1/2) n\ ^ ' 



J {x, y = E = -T^. ^ 122 

^ — ' me *^ — me'2^ 

n—O ^ ' 



This follows from the addition theorem ( 82 §11.4) 



Jy - z) VAtt (t/ + n) r (2t/ -f n) _^ ( \ i ( \ 
(^_,)- = 2-r(. + l/2) 2:- nl ^ Jy+n{w)z J.+„(z) (123) 

Note that the ratio is analytic in w and z for all v. J (x, y) clearly reduces to the previous case 

(f7H|l when = 0. 

It is straightforward to act on this norm kernel with the Hamiltonian to obtain a closed form for the 
Hamiltonian kernel on the dual space, as we did earlier for = 0. We find 

\ J [me " — me^y) 

(me — me^y) 

oo 

= E (i. + nf ^„,, (me-") ^„,, (me^^) (124) 

ri=0 

Once again, this is an hermitian kernel, H (x, y) = H* [y, x), and it smoothly reduces to the flux-free $ = 
case when z/ = 0, as in (|86|) . 

These kernels may be used to compute norms and energy averages, just as in H83|l and H88|) . for arbitrary 
states in the right sector, 

oo 

through the use of the corresponding dual functions 

oo 

V'dual (x) = ("^^") (126) 

Other kernels and corresponding expectation values can be easily constructed, including some like —idyJ (x, y) 
which are non-hermitian. But we will not pursue this here. 
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6 Two exponentials 



This is a complex version of the Morse potential [SS]- With a constant vector potential included, the 
Hamiltonian is 

H — {p + v)^ + /ii exp {ix) + ^2 exp {2ix) (127) 

The analysis of the classical solutions proceeds as before, for fixed energy, and yields trajectories in terms of 
known functions. We forego further discussion of the classical motion (except in a special case given below) 
and consider immediately the quantum system. 

The corresponding Schrodinger energy eigenvalue equation is 



—i— — \- v\ + /ii exp {ix) + /i2 exp {2ix) 



i^E (x) = E^E (x) (128) 



dx 

Single-valued solutions, bounded on the real line, exist for any E > and are given by 

5 WhittakerM f - , ± Ve, 2 J^JI^z ) where z = 

V2 J 

^'--'^^M[^±^-^-^,l±2VE,2V^e''■' ] (129) 




These eigenfunctions are all manifestly even under PT, for fii, and -v/— 7*2 all real. As indicated by the 
notation, there is a double degeneracy for generic E. Here M (a, 6, z) = KummerM (a, 5, z) is Kummer's 
confluent hypergeometric function (pP Chapter 13) also known as i-Fi {a;b;z). 

M (a, 5, .) = (a; b;z) = ± ^^^T^^A (130) 

1 (a 1 (b + n) n\ 

n— 

So the quantum energy spectrum of this Hamiltonian is the positive real line for any choice of /ii and ^2- 
Since it is not always possible to shift x to obtain both real and real /^2, simultaneously, this model is 
not translationally equivalent to a PT symmetric theory, in general, and therefore the model provides an 
interesting counter-example to the naive dictate that real energies require PT symmetry.^** 

The model defines a biorthogonal system for all E' > 0, and for generic E the dual space can be constructed 
using the methods of the previous example involving only one exponential. In particular, if we restrict our 
attention to eigenfunctions periodic on x e [0, 27r], we again find exceptional biorthogonal systems. Due to 
the pre-factor, e(~'^±^/^)"^ 27r-periodicity requires —v± \fE = n G Z. So the energies are given in pairs by 

E±,,^(y±nf for n = 0, 1,--- (131) 

This is exactly the same as in the previous model involving only one exponential. For generic v the 
eigenfunctions again define left and right sectors, as in the single exponential case. The dual space can be 
constructed accordingly. 

If we again restrict our attention to the right sector, for simplicity, the dual space may be constructed 
from polynomials in 1/z. For general [i\ and /i2 the relevant polynomials have been investigated by Erdelyi 
|38j . We defer a discussion of their structure for general [L\ and /i2 until the next section of the paper, where 
we give a unified treatment for all Hamiltonians of the form JQ). For special choices of /^i ^ /i2 , however, 
the biorthogonal structure simplifies dramatically, and many results can again be obtained in closed form. 
If we set = and if we make a particular choice for the coefficients of the two exponentials, somewhat 
surprisingly the eigenfunctions can again be expressed in terms of Bessel functions. 

'^^However, the model is clearly continuously connected to PT invariant theories through independent variations of both ^^\ 
and 112- AH such variations are isospectral deformations of H, i.e. they do not change the total set of energy eigenvalues 
even though they change the form of H and individual energy eigenfunctions. Relatedly, such variations define continuous 
deformations of VT. This statement should be compared to various theorems in the literature on the interplay between general 
anti-linear symmetries and the reality of the energy spectrum |??ll5Hlinnil4| . 
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Simplification for a particular two-exponential case Let 

H = p^ + me" - m^e^" (132) 

Analysis of the classical solutions yields trajectories in terms of elementary functions as in the case of 
one exponential. For _B > 0, the complex momentum trajectories are now "limagons" (for examples, see 
|83|) centered on the Rep axis, rather than the circles of the previous case involving only one exponential. 
There is also considerably more structure in the E = classical solutions. Simple special cases occur when 
g-ix(o) _ ^-^^ jtqj, -v^rjiic}! cases we encounter turning points. 

We now consider the quantum system and immediately note the Darboux factorization |3U| of the 
Hamiltonian. ^ 

H = --^+ me" - m^e^" = ( ^i4- - me' A ( + me'A (133) 
ax^ \ ax J \ ax J 

Schrodinger's eigenvalue equation, the allowed energy eigenvalues, and the periodic eigenfunctions for x S 
[0, 27r] are now given by 

--f^i^n {x) + (me" - m^e^") V'n (x) = n^n {x) (134) 



i^n{x)^Zj^] ^^^^^ 



2fe oo / ,v \ 2fc+l 



^„fc!r(fc + i + n) V 2 ; ,eo^!r(fc 

= Zn\l {/„-i/2 {me") - In+i/2 {me")) 



\/ 2. \ ax J \Jme'^ 

for all integer n, both positive and negative. Z„ is a normalization constant. The /s are modified Bessel 
functions, as given by the explicit series JHSJ, or, because in this case they have semi- integer indices, as given 
by finite sums of powers multiplying hyperbolic functions of me" (jTj §10.2). For n € Z„>o 



The non-degenerate ground state is simply V'n=o {x) — Z^e ™-^" j The result H^n=Q {x) = follows 
immediately from {—ij^ + me") e~""^" = and the factorized form (|133|l . just as in supersymmetric QM 

Duals for all the energy eigenfunctions may be constructed from linear combinations of ■(/'„<o and '0„>o, 
similar to (??). On the other hand, if we restrict to the right sector, i.e. the subspace spanned by 
{V'n I ^ 0}, we may once again make use of only finite polynomials in e~" to construct the dual space. 
Explicit calculation gives 



^ ^ " / (-l)"r(n-fc+l) / me- \ 2''^ ^^"^"^ (-l)^T (n - fc - 1) ^ ^ ^ ^'=+1 ' 

^"("^^-Z,, Ue-j k\ I 2 j + fc! 

= 1 (^TTTT^^-^/^ ^'"^'"^ + (;7^"^"-^V2 {^rne") 

These {%„ | ?i > 0} are biorthonormalized duals for {-!/;„ | n > 0}. By construction, we have 

2 /■27r 



(137) 



^ I Xk (x) ^pj (x) = Skj (138) 



5 The choice Z„ = + n) gives a balanced asymptotic behavior for both bilocal norm kernels, J {x,y) and K{x,y). 
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for all j, fc > 0. The dual functions are now solutions to the inhomogeneous equations 



--^Xn {x) + (me" - TO^e^*'') Xn {x) - ri^Xn {x) = anm^e^" + /3„me*^ 
where the coefHcients of the inhomogeneous terms are given by 



(139) 



an=2k 



(-i)'-+^r(fc + i) 

klZ2k 

(-i)"+ir(fc + i) 



klZ' 



2k+l 



Pn=2k = - {2k+ 1) a„=2fc 
/3n=2fe+l = 2kan=2k+l 



(140) 



The result J^^ Xk (x) ipj (x) dx = for k ^ j again follows from integrating by parts and from J^^ e^^ipj (x) dx - 
— J^^ e^'^Vj (x) dx for all j > 0, as in et seq. 

The construction of the bilocal norm and Hamiltonian kernels proceeds as before. Acting on the dual 
space of the right sector, we find {w = me~", z = me'^) 



r27r 



J {x,y) = X!'^" 

n=0 

d 



(141) 



1 








1 


2^ 






1 


2a 


/wz 



w— w ] \ z— z ] + — 

dw 



dz 



^ |Z„|^ In+l/2 {w) In+1/2 (z) 



7= X] + 1)^+1/2 (w) In+1/2 (z) 



dw 



w + z- z 

dz 



^ |Z„|^ ( n + i ) /„+i/2 (w) /„+i/2 {z) 



ri=0 



The choice |.Z'„|^ = (2n + 1) tt facilitates evaluation of the sums, since ([S2] §11-41 Eqn(6)) 
sinh + z2 — 2wzs Fk I1/2 (Vw^ + — 2wzs) 



\/w'^ -\- z^ 



\/ nP' + z'^ 



2'wzs 



n=0 



— 2wzs V 2 

This resuh, along with (-1)" P„ (s)|,^_i = 1 and 2 (-1)" ^ P„ (s)L=_i 

TT , , ,„x ^ /X sinh (w + z) 
^Y.{n+ 1/2) In+1/2 {w) In+1/2 {z) = ^ 

^ n— 

^ (n + 1/2) n{n + l) In+1/2 (w) In+1/2 (z) = i cosh (z + w) 



J= ^ (n + 1/2) In+1/2 iw) In+1/2 (z) (-1)" Pn is) 

(142) 



-n (n + 1), leads to 



sinh (z + w) \ 2ii;z 



(143) 



{z + w) J {w + zf 

Unfortunately, for the time being, a utilitarian closed form for the last sum in (|141|l eludes us. We have 
only found an integral representation for it (see [72], Vol. II, 5.7.17, Eqn(ll), p 675). 

Given a sufficiently simple result for that third sum, a closed form for the bilocal Hamiltonian H [x,y) 
would then follow by acting with H on J [x, y), as in the previous cases involving one exponential. 

Canonical integral representations for the wave functions, and their duals, can also be exhibited for the 
two exponential model, especially in the simplified case (|132ll . The easiest of these is probably Poisson's 
integral (see ^ 10.1.13 and 10.1.14 with z iz), or alternatively just H47(l with ^/E replaced by n + 1/2 
and m —f im in the argument of the Bessel function. Corresponding integral representations for the duals 
follow from biorthonormality, for given choices of the kernels (e.g. for Poisson's integral, an appropriate 
kernel is given by 10.2.36). This allows construction of a Table as in the case of the one exponential 
model. We leave the details of this construction to the well-motivated, interested reader. 
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7 Biorthogonal systems for H = + X^fc>o f^k exp (ikx) 

We indicate methods to analyze the general case, reserving some of the details to be presented elsewhere. 
For X g (—00, +00) without any periodicity constraints, the spectrum of this Hamiltonian is known j42) 
to be the real, positive half-line, for any choice of /zs such that J2k>o < 00, and not just for those /zs 
which are real or a;-translationally equivalent to real values. So PT symmetry is not required for real energy 
eigenvalues in this model. Here, we will restrict our attention to 27r-periodic functions and their duals, to 
obtain a discrete subset of real energy eigenvalues, namely just {i?„ = | n = 0, 1, • • • }. 

A simple procedure to construct this biorthogonal system is to begin with the dual functions. We 
assume these to be polynomials in powers of = e~", so we write 

Xn (-) = 37 E ' H=(z-] l^uz'^ (144) 



j=o ^ ' fe>0 

These polynomials, for three or more /is, would seem to be a natural generalization of those constructed by 
Neumann, Gegenbauer, and Erdelyi. Acting with H — v? gives 

^ fc>0 i=i fe>o j=Q 

We now impose the condition that the RHS involve only positive powers of z so as to obtain an inhomogeneity 
that will be orthogonal to the span of {z" | > 0} under contour integration • • • , as in previous examples. 
This leads to n equations that fix the coefficients Cn,k for A; = 1, • • • , n in terms of c„fl, the latter being an 
overall choice of normalization. Explicitly, the powers to be eliminated, and the resulting "triangular" set 
of n equations, are 

: ({1 - nf - nA + mcnfi = 

Z""+2 : Cn,2 ( (2 - n)^ - ) + /Lt2C„,o + MlCn,l = 

^ ^ (146) 

^-n+n . p^^^ ^(^^ „ _ ^ ^InCn.Q + Mn-lC„4 H h /^lC„^„_i = 

These equations can always be solved, sequentially, for c„.k in terms of Cn,o and the /is, as expressed by the 
recursion 

1 ''"^ 

"""■^ " (2n-k]k ^ ^^-■'■^"^J ' fc = 1, • • • , «• (147) 

So for a given n the Cn,k will depend on /ij<„ but not on ^j-^n- Indeed, the c„^fc are finite polynomials in 
the /iS, and therefore non-singular for all /is. 

Thus, up to their individual normalizations, the dual functions are completely determined along with 
the inhomogeneous equations that they obey, without (yet!) having to find the eigenfunctions of H. The 
resulting equation giving the action of H on the dual functions is of the form 

(i7-n2)x„(z) = E7n,fe^' (148) 

k>Q 

where the coefficients ^n,k may also be expressed explicitly as finite polynomials in the //s, 

n 

ln,k ^^fJ.k+7i~jCn,j (149) 

j=0 

^^So these systems are an exception to the statement: "A direct and frontal attack on the problem of determining [the dual 
polynomials] is not usually fruitful." - Morse and Feshbach I56l p 931. 
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This result follows from H145|l and (|147|) . Note that k is not necessarily limited by n in 7„.fc. The 7„^fc 
will depend on all the /is, in general, with k taking on all values up to and including the highest power of z 
appearing in H. 

With {xn (z)} in hand, we may then proceed to construct {ipn (z)} as analytic functions around the origin 
by demanding biorthonormality with {x„ (z)}. We take the ipn functions to be infinite series of positive 
powers 

oo 

If we assume the requisite convergence of these series, their forms automatically give §^Xk {z)ipn {z) = 
for n > k. Otherwise, conditions are obtained on the coefficients a„j for j + n < k hy requiring 
^Xk (z) Ipn (z) = Skn for cvery k > n. These conditions amount to another triangular set of equations 
which can always be solved, sequentially, for the a„j-, j < k — n, for any choice of ^s. Namely 



1 

Cn.O 



k~r 

E 

3=0 



Ck,k- 



■n—j ^n.j 



for k > n 



(151) 



By considering all fc > n in succession, we obtain all a„j . The series for ipn is a development in the minors 
that invert these triangular equations. For convenience we choose Cn,o = 1 for all n, then 



Xn (z) = 4r (l + Cri,lZ + C„.22^ 



iPn (z) = Z" - C„+1,1Z"+1 + 



Cri+1,1 1 
Cn+2,2 C„+2,l 



„n+2 



Cn+ia 1 
Cn+2,2 C„+2,l 1 
Cn+3,3 Cn+3,2 C„+34 



(152) 
(153) 



The coefficients of in -pn (z) are again finite polynomials in the /iS. While convergence of this series 
is certainly not obvious for arbitrary /is, as written, it is clear that convergence can be determined on a 
case-by-case basis from the explicit form of the coefficients. Moreover, when the number of /is is finite, no 
matter what their values are, it is not too difficult to show the ipn {z) are entire functions of z. 

Thus, if we assume the requisite convergence of the ipn series, all the {ipn {z) | for n > 0} are determined, 
and they are complete on the span of {z" | n > 0}. All positive powers of z can be expressed as series of 
{ipn {z)}, similar in form to H20|l . just as all negative powers of z can be expressed as finite sums of {xn {z)}, 
similar in form to l|21(l . 

Remarkably, the {ipn {z){ just obtained turn out to be non-degenerate energy eigenfunctions, and are in 
fact all of the 27r-periodic eigenfunctions of H , with 



Hipn {z) = n^lpn {z) 



(154) 



The RHS of 1)154(1 can be deduced in a novel way from p48|l . the biorthogonality of {xj {z) ,ipk (z)}, and 
the completeness of {ipn (z)} for analytic functions about the origin, as follows. Given that ipn {z) is such 
an analytic function, as in (|150l) . completeness allows us to write Hipn (z) = X^fcLn ^n,k'>pk (z). From this 
and biorthonormality, we have bn^k — '^§^Xk {z) Hipn (z). But then, upon integrating by parts and using 
(|148|l as well as orthogonality relations, we also have -^§^Xk {z) Hipn {z) = n'^Sk,n- So bn^k = n'^Sk,n, and 
(|154|l is obtained. ^"^ 

Conversely, given ((154|l and 1(148(1 . we may prove the orthogonality relations §^Xk {z)ipn {z) = for 
k ^ n just by inserting H and following the same line of argument as in the case of a single exponential 
potential, 1(54(1 et seq. Some explicit closed forms for the various functions, some kernels constructed from 
sums of their bilinears, and related convergence issues, as well as many other details for this general case, 
will be discussed elsewhere. 



^^Note that 2 = is a regular singular point of 11541 . for any number and any values of the ^s such that X]fe>o Mfe-^*^ converges 
near the origin. In fact, 11531 is exactly the conventional series obtained by expanding about the regular singular point at 
z = 0, albeit the series was obtained here in an unusual way from the properties of the dual space polynomials. Convergence 
criteria for 11531 in terms of the /js are obtained in the conventional series approach in 1731 1421 . without any discussion of the 
dual polynomials. 
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8 A brief look at QFT extensions 



For the complex Liouvillc quantum field theory in 1+1 spacetime, we choose to define the Hamiltonian 
density by normal-ordering at mass m in the standard way [23]. 

n = N„, Qtt^ + i {d^^f + f,e^^P^^ (155) 

Normal-ordering at any other mass should not affect our conclusions. The resulting (pseudo)scalar field 
theory is conformally covariant [291 , at least for a range of /?. For this field theory, the only issue to 
be raised here is the stability of the ground state. HoUowood has previously discussed separate but 
related issues in his study of solitons for complex Toda theories using conventional Hilbert space methods 
and perturbation theory. 

We may attempt to place an upper bound on the ground state (7i) using the conventional Hilbert space 
norm through an adaptation [23 of the variational argument of §111 in Coleman [53| . Take as a trial state 
the zero- mode-shifted Fock vacuum for arbitrary mass M , defined by the effect of annihilation operators for 
that mass: 

a(0,M) |0,C,M) = i,^|0,^,M) , a (fc, Af) |0, ^, M) = if k^O (156) 

For |0, ^, M) the zero-mode configuration is coherent, while the a{k,m M) non-zero mode configuration is 
superfluidic (i.e. built formally from powers of correlated pairs of creation operators {a) {k,m)a^ {—k^mj) 
acting on |0,0, m)). Note that 

(0, ^, M| iV™ {e"'^^) |0, C, M) = e2./3C (q, M\ N„, (cos (2/3 (V- - 0)) |0, ^, M) (157) 

since (0, M\ (sin (2/3 (tjj - C))) |0, M) = 0. For this trial state, we find 

(0,e,M|H|0,C,M)=/.e2»^« — + (M'-m^) (158) 

\ J Sn 

Therefore, by choosing ^ such that /xe^''^^ is real and negative, and then letting M/m oo, there is no lower 
bound on (0, ^, M\ TL |0, ^, M) if /3^ > 27r. Thus in the usual Fock space framework, there is no ground state 
if 0^ exceeds 27r. This may indicate a phase-transition for the field theory. 

For the field theoretic extension corresponding to the model with v — the RHS of (|158|l becomes 
just a sum over k ~ 2(3 with various /i^.. Thus the variational argument shows an instability whenever the 
sum over k includes fc > 5. 

The problem is now to determine how this variational argument is effected by using the bilocal norm 
instead of the usual Hilbert space norm. The operator methods in 118' and/or the functional methods in 
|27j should suffice for this purpose. There have been studies of the variational principle for similar "quasi- 
hermitian" theories [J^ and the techniques introduced in those studies should be applicable to expectations 
of (|155|l . This issue will be fully discussed elsewhere. 

Here we only remark that the situation is reminiscent of —(j>* field theory, which also has no ground state 
based on conventional variational arguments using the standard Hilbert space norm and trial states with 
real expectations {(f)), but which appears to be an acceptable model when viewed as PT symmetric cjP {icff 
theory ^Uj. This is unquestionably related to the use of the CPT norm, but more directly, it is due to the 
definition of the zero- mode problem on a contour in the LHP, and not the real axis. This obviates most 
physical intuition about the ground state of that theory, in our opinion. Even so, there are obviously major 
differences between PT symmetric — </>^ field theory and imaginary Liouville theory as defined by l|155l) . in 
that the apparent instability in — (/)^ field theory is evident at the classical level for any coupling strength, 
whereas the instability indicated by H158|l for 0^ > 27r is a purely quantum effect. 

9 Conclusions 

We made a thorough analysis of some complex one dimensional periodic potentials, of Liouville type, within 
the framework of biorthogonal systems. We provided an extensive guide to the considerable literature for 
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these models. We conclude that these systems are consistent at both classical and quantum mechanical 
levels, with remarkable and deep connections between the two levels, just as in the case of real Liouville 
models. In fact, the connections are more intimate in the complex case, inasmuch as the classical structure 
is richer. We studied in some detail the correspondence and impact of the classical dynamics on that of the 
quantum problem. 

More specifically we have shown that the complex Liouville model is equivalent to a projected free 
particle system, and we may choose that free particle projection to be chiral. With that choice, the complex 
Hamiltonian may be viewed as a completely local description for the free chiral particle, the only such local 
Hamiltonian description known to us. 

We have used the remarkable properties of Neumann polynomials as the biorthogonal partners of Bessel 
functions on a circle in the complex plane. With those properties in hand, we obtained closed forms for 
the inner product and various other hermitian kernels. We enlarged this study to include a wide class 
of complex periodic one dimensional Hamiltonians, and we obtained similar closed form results in several 
additional cases. 

In our analysis, the theory of PT symmetric models was often a guide but did not play a decisive 
role insofar as the states were self-orthogonal in the periodic sectors that we emphasized and some of the 
models that we considered were not PT symmetric. Moreover, many of the concepts that we used are already 
contained in the previously existing framework of biorthogonal systems. As stressed in the title of the paper, 
the use of biorthonormal systems of vectors for the construction of consistent quantum mechanical systems 
was the real conceptual lift, applicable even when the states are self-orthogonal. The most important feature 
shared with the PT symmetric approach was simply that of a complex Hamiltonian with real eigenvalues. 

It would be interesting to investigate biorthogonal systems for self-orthogonal subsectors of other solvable 
complex potential models, e.g. those in for which the duals would seem to be the associated functions 
described in [HOI- Another interesting and possibly important challenge is to combine the biorthogonal 
approach with the study of symmetries of the Hamiltonian. Beyond the usual commutators between any 
local symmetry generators, there must also be commutation relations among the various bilocal kernels. 
Thus the full set of symmetries of the corresponding quantum problem may not be local and may not be 
obvious. Along these lines, it is perhaps a useful working hypothesis that any consistent local complex 
Hamiltonian with real eigenvalues is related by a similarity transformation to that of a hermitian, albeit 
nonlocal, Hamiltonian. However, this conceptual device may not always lead to practical calculational tools 
as this similarity transformation is usually not known, a priori, and may be difficult to implement when it 
is known. 

A clear and readable discussion on the use of similarity transformations to expose the hermitian un- 
derpinnings of certain apparently non-hermitian (i.e. "quasi-hermitian" ) theories is given in }75J (also see 
|(i3l lf)4l Ififil \(>7\ ). The general theory is nicely illustrated in a recent paper on PT-symmetric —z'^ quan- 
tum mechanics ^Hl; where after the change of variable z = —2i\/l + ix, x £ M., the relevant similarity 
transformation is very simple and easily implemented. 

While substantial progress has already been made on several PT symmetric field theories [S] |H7], it 
is important to apply the methods of biorthogonal systems to more general quantum field theories and 
many-body problems, particularly those of the complex Liouville and Toda type. It remains a challenge to 
consider situations involving self-orthogonal states in these models, to compute correlation functions using 
those states, and to investigate their possible relevance for string theories. 

The study of the present paper may lead to a deeper appreciation of the behavior of quantum Liouville 
field theory for complex coupling constants, as sketched in Section 8. This would further understanding of 
subcritical string theory, especially the correlators for those models which exhibit interesting structure in the 
complex coupling plane (cf. [021 and references therein, especially (HEIEZI)- In fact, Liouville theory with 
imaginary coupling (or its close cousin under complex continuation - so-called "timelike" Liouville theory) 
has already been analyzed using conformal field theory methods [7Hl[7n|. We believe the Hamiltonian analysis 
initiated above will help shed additional light on this subject. 

Moreover, the field theory extensions of the models described in Section 7 are relevant to applications of 
conformal field theory and phase-transitions. In particular, complex Liouville field theory may be useful to 
understand freezing transitions j21l I22| . More recently, there is very interesting work 2 involving complex 
Liouville field theory in non-Hcrmitian Luttinger liquids and flux line pinning in planar superconductors. As 
for other possible physical applications of the techniques in this paper, it is perhaps interesting to conjecture 
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that coupled Bose-Einstein systems with dissipative effects might be governed by relatively simple non- 
hermitian effective Hamiltonians. This could be, in several respects, a completely quantum analogue of the 
classical driven/damped oscillator systems discussed in |I7) . 
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